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Abstract 



The contributions of virtual corrections and soft gluon emission to the inclusive Higgs boson 
CNj ' production cross section pp H + X are computed at next-to-next-to-leading order in the heavy 

top quark limit. We show that this part of the total cross section is well behaved in the sense of 
, perturbative convergence, with the NNLO corrections amounting to an enhancement of the NLO 

' cross section by ~ 5% for the CERN LHC and 10-20% for the Fermilab Tevatron. We compare our 

■ results with an existing estimate of the full NNLO effects and argue that an analytic evaluation of 

the hard scattering contributions is needed. 

o 



Cli! 1 Introduction 

(— I I In the past two decades the predictions of the standard model have been confirmed with 

remarkable precision. Still, the agent of electroweak symmetry breaking remains elusive. 
. The simplest mechanism introduces a single complex scalar doublet. Three components of 

r> I the doublet give mass to the W and Z bosons, leaving a single neutral particle, the Higgs 

' boson, as the signature of the symmetry breaking sector. This scenario is called the min- 

imal standard model (SM) and forms the benchmark for the investigation of electroweak 
symmetry breaking. 

In their final run, the experiments at the CERN e+e~ collider LEP established a lower 
mass limit for a SM Higgs boson of ^ 113 GeV [Q]. This value severely challenges the reach 
of the Fermilab Tevatron. Nonetheless, there are indications that the Higgs boson is not 
much heavier than this limit. Fits to precision electroweak data actually prefer a value of 
the Higgs boson mass that is well below the exclusion limit and place a 95% confidence 
level upper limit near 200 GeV Q. Supersymmetric extensions to the SM, which have at 
least two Higgs doublets and therefore at least five observable Higgs bosons, prefer that the 
mass of the lightest Higgs be below approximately 135 GeV |^. Detailed studies indicate 
that, with sufficient luminosity, the Tevatron could be sensitive to Higgs bosons somewhat 
beyond the W^W~ threshold [§. 
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At hadron colliders, the dominant production mechanism for Higgs bosons of mass below 
~ 700 GeV is gluon fusion. The Tevatron however will only be able to make limited use 
of this production channel. For Higgs boson masses well above the W^W~ threshold, the 
production rate is too small to be observed. For masses well below the W^W~ threshold, 
the dominant bb decay mode is overwhelmed by QCD background, and the production 
rate is insufficient to allow observation of rare decay modes like ^ 77. Thus except 
for a window around the W^W~ threshold where that decay mode may be observable, 
the Tevatron Higgs search will rely on associated production with aW 01 Z boson rather 
than gluon fusion 

At the CERN Large Hadron Collider (LHC) the gluon fusion production mechanism will 
be extremely important. For lighter Higgs boson masses, the bb decay mode will still 
be overwhelmed by QCD effects, but because of the high machine luminosity and the 
large gluon luminosity at relatively small parton energy fractions, the LHC will be able 
to measure the ^ 77 decay mode. This will be important not only as a discovery 
channel (if the Higgs boson has not yet been found) but also as a means of studying Higgs 
properties like its coupling to top quarks. 

Unfortunately, the gluon fusion process is currently not fully under control. The next- 
to-leading order (NLO) corrections to the production rate were computed ten years ago 
and were found to be extremely large, of order 70 — 100% Such large corrections 

clearly ask for the evaluation of higher order terms in the perturbative series in order to 
arrive at a solid theoretical understanding of the process. In this paper, we present the 
soft plus virtual corrections to inclusive Higgs production at next-to-next-to-leading order 
(NNLO). For the relevant values of Mh ~ 100-200 GeV and the center-of-mass (c.m.s.) 
energies at Tevatron (2 GeV) and LHC (14 TeV), these terms are not expected to dominate 
the full result. Nevertheless, they represent a first step towards the complete answer, in 
the sense that they comprise a well-defined, ultraviolet and infrared finite, gauge and 
renormalization group invariant piece. By comparing to existing estimates on the NNLO 
corrections for gg ^ H, we conclude that the evaluation of the full answer is necessary. 

2 Higgs production through Ught parton scattering 
2.1 Leading order 

Assuming the standard model as the theory of particle interactions, the coupling of the 
Higgs boson to gluons is mediated by a quark loop. In the following we will neglect all 
quark masses except for the top mass. In this case only the top quark contributes, because 
the Yukawa couplings are proportional to the fermion masses. The lowest order diagram is 
then shown in Fig.||(a). It was computed some time ago and leads to the cross section 
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Figure 1: Leading order diagram to the process gg H: (a) in full QCD; (b) in 
the effective theory [Eq. (Q)]. The straight solid lines represent the top quark, the 
symbol ® denotes the effective vertex. 



where s is the partonic center-of-mass energy, G-p the Fermi coupling constant, Us the 
strong coupling constant (which depends on the renormalization scale ^), Mt the pole 
mass of the top quark, and Mh the Higgs boson mass. 



2.2 Effective Lagrangian 

The lowest threshold of the diagram in Fig. |^(a) is at Mh = 2Mt. This means that an 
expansion in terms of M^/Mt is expected to converge for Mh < 2Mt. In fact, it has been 
shown that the NLO K-factor for this process is excellently approximated up to much 
larger Higgs boson masses, even by keeping only the first term of such an expansion. With 
a top mass of around Mt = 175 GeV and the experimental data favoring a relatively low 
Higgs boson mass between 100 to 200 GeV, we feel that keeping only the leading term in 
Mn/Mt is also well justified at NNLO. 

We therefore integrate out the top quark and compute amplitudes using QCD with five 
active flavors and the following effective Lagrangian for the Higgs-gluon interaction: 

= Oi = -£c7f Of, Oi = G^G'^'^'' , (2) 

where is the gluon field strength tensor. In the approximation that all light flavors are 
massless, this effective Lagrangian is renormalization group invariant, but the coefficient 
function Cf and the operator Of must each be renormalized. Below, we give the renor- 
malized value of Ci while the renormalization of Of and its matrix elements is discussed 



in Sect. 3.1 



The coefficient function contains the residual logarithmic dependence on the top quark 
mass and has been computed up to 0(af) |l^. For our purposes, we need it only up to 
0{Q^g) ||lC,ll]. The coefficient function in the modified minimal subtraction scheme (MS) 
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(3) 



where = ln(^^/M^^) and Mt is the on-shell top quark mass, as = ai^\fi'^) is the 
MS renormalized QCD coupling constant for five active flavors, and ^^/ is the number of 
massless flavors. In our numerical results, we always set nj = 5. 

Using the effective Lagrangian of Eq. (^) , the lowest order amplitude for the process gg — > 
H reduces from the one- loop diagram in Fig.|l](a) to the tree-level diagram in Fig.[l|(b). 
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Figure 2: Sample two-loop diagrams contributing to the virtual corrections at NNLO. 
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Figure 3: Typical diagrams contributing to (a) single and (b) double real radiation of 
gluons and (c) double real radiation of quark-antiquark pairs at NNLO. The dashed 
line is the Higgs boson and "(8)" denotes the effective Higgs-gluon vertex. 



2.3 Sub-Processes: virtual, single real, double real 

For a complete treatment of the process pp ^ H in NNLO QCD one needs to compute the 
partonic processes which must then be convoluted with parton distribution functions. 

The partonic subprocesses are: 

(i) Virtual corrections to two loops: 
99 ^H. 

(m) Single real radiation to one loop: 

99 ^ Hg, qq^ Hg, qg ^ Hq, qg ^ Hq. 

{in) Double real radiation at the tree level: 

99 Hgg , gg Hqq , qq Hgg , qq Hqq, qg Hqg , qg Hqg. 

2.4 The soft approximation 

Radiative corrections to inclusive production fall into three categories: delta function 
terms cx (5(1 — x), large logarithms of the form ln"(l — x)/(l — x), and hard scattering 
terms that have at most a logarithmic singularity as x — > 1, where x is the ratio of the 
Higgs boson mass squared to the incoming parton c.m.s. energy squared [see Eq. (|^)]. 
In the kinematically similar Drell-Yan process and in inclusive Higgs production at 
NLO 1^], it has been observed that the delta function and large log terms are important 
though not necessarily dominant parts of the radiative corrections. 

Virtual corrections add only delta function contributions while real emission corrections 
contribute to all three categories. The soft approximation is obtained by taking the limit 
X ^ 1. In doing so, one is left with only those terms that scale like (1 — x)"^""*^, where 
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e = (4 — (i)/2, and d is the number of space-time dimensions assumed for the regularization 
of divergent integrals. The term (1 — x)"^""*^ must be interpreted as a distribution, 
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where Vnix) is a "plus" distribution, 
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(5) 



Thus, towers of the large logarithm terms 'Dn{x) can be trivially obtained from knowledge 
of the 5(1 — x) terms. 



3 Results 

3.1 Renormalization of matrix elements 

We write the partonic cross section as follows: 

vr 1 
^ = ^^^1 («s) _ + PB + Pc] ■ (6) 

PA is the contribution from virtual corrections, while pB and pc represent the contributions 
from single and double real emission, respectively. Z\{a^ is the global renormalization 
factor for the composite operator in the effective Lagrangian p^: 



2 fl-n „ (7) 

+ 0(a3). 



a2 



with Zi(q^ = — — — = 1 + — 

1 — p(as)/e vr e V vr / 

Here, /3(as) is the QCD beta function: 

/3(a,) = -^/3o-(-)'/3i + 0(a3), 

1 / 2 \ n ^ ( 38 \ 

/?o = ^(ll--n,j, /5, = -(l02--n, 

In the following sections, we will list the results for p\^ ps, and pc in terms of the bare 
coupling which is related to the MS renormalized coupling by 

where = exp[e(— 7e + ln4vr)] , 

and 7e = 0.577216 ... is the Euler constant. The charge renormalization constant Za{as) 
is 

^«(a.) = l- --+(-) (10) 
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3.2 Two-loop virtual corrections 



The two- loop virtual corrections have been calculated in [|T^ . Typical diagrams are shown 
in Fig. ^. Since gluons are massless, these diagrams depend on only a single scale = Mfj, 
where q is the momentum of the Higgs boson. They were calculated using the method of 
ref. [f^ which allows one to relate the integration-by-parts identities for the two-loop three- 



point functions of Fig. g to those for three- loop two-point functions |16|. The coding of the 
integration routines was done in FORM |17| and was based on the program MINCER [18|. 
It is worth mentioning that the resulting program is capable of dealing with much more 
general diagrams than the ones considered for this calculation. In particular it solves 
integrals with (in principle) arbitrary integer powers of the propagators. This may become 
useful in future applications where expansions have to be applied to the diagrams. 



The expression for is given in ref. 
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where Cn = C("') is Riemann's C function which takes the particular values 

(2 = 7rV6 = 1.64493... 
Cs = 1.20206... . 



(12) 



(13) 



3.3 Single real radiation 

The one-loop amplitudes for the processes {ii) in Sect.^]^ were computed in ref. |]l9| 
in the form of helicity amplitudes to order e''. Unfortunately, these amplitudes are not 
sufficient for our purpose. The integration over two-particle phase space generates infrared 
singularities which take the form of single and double poles in e which combine with 
the 0{e) and O(e^) pieces of the amplitude to produce finite contributions to the cross 
section. If we were to work in the four dimensional helicity scheme [^], where gluons have 
two polarization states, helicity amplitudes extended to higher order in e would suffice. 
However, since we work in the conventional dimensional regularization scheme, gluons 
have d — 2 = 2 — 2e polarizations. Thus, in addition to helicity amplitudes we also need 
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amplitudes corresponding to the extra "e" polarizations. Both the additional pieces of 
the amplitudes and the integration over the two-particle phase space have been evaluated 
and will be presented elsewhere [^] . That calculation includes the full dependence on the 
center-of-mass energy. Here we are interested in the soft limit, to which we find that only 
the gg Hg process contributes. The result is 

PB = fp^^ + [f)\^^+Oial), (14) 
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Expanding fae{x) (see Eq. (^), Eq. (|15|) describes both delta function and large logarithm 
contributions to the radiative corrections. 



3.4 Double real radiation 



The tree-level amplitudes for the processes {in) of Sect. 2.3 with two partons in the final 
state was evaluated in |2^, though again only as helicity amplitudes to 0{e^). We have 
evaluated the full e dependence in the conventional dimensional regularization scheme by 
computing all of the Feynman diagrams, squaring the amplitude and integrating over phase 
space, dropping terms that do not contribute in the soft limit as described in Sect. p.4| . 
We find that only the gg Hgg and gg — > Hqq processes contribute with the result 
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The terms proportional to n/ come from the gg Hqq sub-process. (See Fig.|3[ 



3.5 Mass factorization 



After combining all processes and renormalizing, there are still infrared singularities left 
over. These are all associated with mass factorization and can be absorbed into process in- 
dependent functions Tij{x) associated with the incoming partons. The (IR-finite) partonic 
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cross section, a, is defined implicitly by the relation 

CTij = r 



(17) 



where a is defined in Eq. (P). The symbol (8" indicates convolution over longitudinal mo- 
mentum fractions, 



if^g) (x) 







dy / dz f{y)g{z)6{x -yz), 



(18) 







and the subscripts z, i, etc., indicate parton identities. The Tij{x)^s are given by 



Tij{x) = 5ij6{l - X 
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where the p/j^'' are the nth order splitting functions |23,24]. In the soft limit, only pj^J^^ 



with i = j = g contribute. They are given by 
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and similarly for agg, Eq. (|T^) can be solved for agg order by order to give 
^!'J{x) = afJ{x)^ao6il-x), 
aS{x)=a('J{x) + 2aoP!j^\x)-^, 

^SH^) = - (/5o + 2 (Pj? Pi°))(x 

+ ^(aoP«(x) + 2(p(o)®a«)(x) 

where 
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The convolutions in Eq. (^) can be evaluated with the help of the following relations 
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C|5(1-X)-2C2P1(X)+P3(X). 



3.6 Partonic cross section 



Putting things together, we finally arrive at 
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with Ijj = /Mfj) and = In(;u2/M(). As a numerical example, let us set /x^ = Mfj 
and insert the values Mh = 130 GeV, Mt = 175 GeV, and a., = 0.112: 



vr /Q 



0", 



99 



576^2 \ vr 



[1 + 0.548 + 0.107] 5(1 - x) + [0 + + 0.283] Vq{x) 
+ [0 + 0.428 - 0.040] Vi (x) + [0 + - 0.029] V2 (x) (26) 
+ [0 + + 0.092]p3(x)l +0(af), 



where the individual numbers in square brackets correspond to the d'^gg, o'^gg, and afg 
contributions. We observe that the coefficient of the 5 function exhibits satisfactory con- 
vergence behavior. The magnitude of its NNLO contribution agrees roughly with the 



estimate of ref. [14|. 



3.6.1 Checks on the result 

There are several checks that can be performed on our result. The first check is to observe 
that all poles in e cancel to give us a finite result for a'f'g'^"^^ . Since the Vn{x) terms at 
order e'^ are linked to the poles of the real emission contributions, this already gives us 
great confidence in these terms. Below we will discuss an even more stringent constraint 
on the T>n coefficients. 

Another check is to compare the three leading poles of the virtual two-loop corrections 



to the general result of ref. |25|, where they are expressed in terms of universal functions 
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that depend only on the identity of the external partons. The poles of the diagrammatic 



result of 1 14] fully obey this observation p6| . 

One can now combine the universal (three leading) NNLO pole terms with the NLO results 
yO^'^^ and to derive the coefficients of Vn{x) for n = 1,2,3 in a^gg'^°^^ . Note that 
this uses only one-loop results and the universal behavior of QCD amplitudes and is done 
without reference to any of the newly calculated NNLO terms (including those of p^). 
The pole structure of the real radiation terms must take the form 
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Expanding the fne{x) terms as in Eq. (^ and requiring that all poles vanish from al^J , we 

can solve for an4, fflns and an2 (n = 2,4), while 021 and 041 remain undetermined. This 

(2) 

solution is sufficient to fix the coefficients of I'i(x), 1^2(2;) and 1^3 (x) in agg , though not the 
coefficient of T>o{x). We find complete agreement between the T>^^ 2,i} terms determined 
in this way and those found in Eq. (p5|). 

Let us note that the absence of terms oc fne with n 7^ 2, 4 in Eq. ( p7D can be understood 
from the following consideration. The index n of the fne terms has two sources: the phase 
space elements for single and double real radiation contribute factors of (1 — j;)^"^*^ and 
(1 — x)'^"^'^ respectively (see e.g. Appendix E of ref. fl^), while the one-loop corrections 
to single real radiation amplitudes squared can contribute factors of (1 — x)"*, (1 — x)™"*^ 



or (1 — x)"^ , where m is an integer [19, 21|. [Tree-level amplitudes squared contain 



only integer powers of (1 — x).] The (1 — x)~ terms in the one-loop amplitude squared 
should generate an /3e(x) contribution. However, the universal functions describing the 
soft factorization properties of one-loop amplitudes only contain terms that scale like 
(1 — x)~^ and (1 — x)"^"^*^ [^. Thus, there are no terms in the squared amplitude that 
scale like (1 — x)~^~'^ and that survive the soft limit. 

Another way to derive the 'D{i^2,3}i^) terms without using any process-specific NNLO 
results is to perform a threshold resummation analysis along the lines of [^,|ll[. Again 
we find complete agreement with our result. As a by-product, we can use our result for 
the coefficient of Pq in order to derive additional input for the resummation formula. Of 
particular interest is the NNLO coefficient attributed to large-angle soft-gluon emissions 
which has been derived for the Drell-Yan and deep inelastic scattering processes in ref. [29 
(where it is called -Df , P = DY,DIS, . . . ). Following this analysis for gg — > H, one finds 
that the corresponding coefficient can be obtained from the one in the Drell-Yan process 
[Eq. (21) in |2^] simply by changing the color factor Cp to C^-Q 

A further check concerns the terms proportional to In //^ which can be derived from renor- 
malization group and factorization scale invariance of the physical cross section:^ 

CTpp = agg g ^ g , (28) 



^We thank W. Vogelsang for pointing this out to us. 
^Note that we restrict ourselves to the limit of soft gluons. 
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where g is the gluon distribution function. This leads to 



Equation (p5|) satisfies this relation. 



^ 99 ' \ ^ I 99 



agg = 0{al) . (29) 



4 Discussion 

For a consistent discussion of the physical implications at NNLO, one has to convolute 
the partonic cross section at a given order with the corresponding parton distribution 
functions (PDFs) of the same order. Unfortunately, NNLO PDFs are not yet available.^ 
Thus, we will restrict ourselves to a semi-quantitative investigation of the NNLO effects 
by using NLO PDFs at NNLO. In ah results we have used the CTEQ5 ^ family of parton 
distributions. 

We define 

-f^NLO = ^^"^ , -fi^NNLO = ^'^^'^ , (30) 
O'LO ctlo 

(2) 

where (Tnnlo is obtained by convoluting the partonic agg with NLO PDFs and using NLO 
evolution of the strong coupling constant, ctlo and o"nlo are evaluated consistently with 
LO and NLO PDFs and evolution, respectively. Of course the treatment for i^NNLO is 
inconsistent, but one may nevertheless get an idea on the magnitude of the effects. 

Figure |^ shows the i('-factors Knlo and Knnlo for a proton-proton c.m.s. energy of 
(a) 14 TeV and (b) 2 TeV, when only the purely soft terms [i.e. oc (5(1 — x) and oc Vni^x)] 
are kept. One finds a nicely convergent behavior when going from NLO to NNLO, both 
for LHC and Tevatron energies. Also shown are the approximate soft NNLO terms of [|ll| 
(see below). 

Figure |5| shows the cross section a{pp H + X) obtained by weighting the terms in 
Eq. ( p5|) by the ratio of the LO expression for a including full top mass dependence given 
in Eq. (|^) with the LO result in the heavy top limit. This procedure has been shown to 
be an excellent approximation of the full top mass dependence at NLO for Higgs boson 
masses up to and even beyond the top threshold. These curves again demonstrate the nice 
convergence of the soft contributions to the cross section. 

Note that these numbers do not represent the full phenomenological result, because of the 
missing, non-negligible contributions coming from hard scattering. Rather, this consistent 
comparison of the soft limits at NLO and NNLO should give an indication of the magnitude 
of the total corrections. 



4.1 Estimating the hard scattering contributions 

At NLO one finds that the soft corrections, while significant, are by no means dominant. 
Comparison of the dashed lines in Fig. ^ to those of Fig. ^ (discussed below), shows that 

■^Various necessary ingredients are already available, and it seems to be clear that the combined efforts 
of different groups will eventually result in an appropriate set of distributions. 
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Figure 4: K factor as defined in Eq. (|^), using the purely soft approximation of 
Eq. (25). Dashed and sohd hnes correspond to NLO and NNLO, respectively. The 



dotted line represents the approximate result Ug"^* of |11|. 
(a): ^^S= 14 TeV, (b): ^^S = 2TeV, where ^^S is the c.m.s. energy of the proton- 
proton system. 
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Figure 5: Cross section a{pp H + X) in the purely soft approximation (c/. 
caption of Fig.|^. Dash-dotted, dashed, and solid line correspond to LO, NLO, and 
NNLO results obtained from Eq. (|25|). The dotted line is the soft part of the NNLO 
approximation of ref. Q. (a): VS = 14 TeV, (b): VS = 2 TeV. 
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the soft terms account for no more than half (and often significantly less) of the total NLO 
correction. Thus, it will be necessary to compute the full NNLO correction, including the 
hard scattering terms, before one could conclude that the gluon fusion process is under 
control. 

Until such time as the full calculation is available, we rely on the estimate of ref. |]ll| 
which points out that the hard scattering corrections at NLO are dominated by formally 
sub- leading terms of the form ln*(l — x) and then estimates the NNLO correction including 
these sub-leading terms. We write their result as 

Ggg is generated from an expansion of a one-loop resummation formula and naturally the 

iKrf in fnrm of ?, 

99 



analytic form of 0"^°^* differs from the actual NNLO result given in Eq. (25). Similarly, the 



(formally sub-leading, but numerically dominant) part is only an approximation of the 
actual, as yet unknown expression. It is expected that the terms with the highest power 
of ln(l — x) [ln(l — x) at NLO, ln^(l — x) at NNLO] obtained in this manner are correct but 
that terms with lower powers of ln(l — x) are not. In the case of Drell-Yan production. 



the result of ref. |11| including the corresponding sub-leading terms cr^g reproduces the 



full answer to within less than 5%. However, as can be seen in Fig.^, the soft terms 
generated for Higgs productions do not agree with the true NNLO soft terms as precisely, 
especially at LHC energies. 

With these caveats, we present an estimate of the full result by combining our result for the 
soft correction with the approximate, formally sub-leading logarithmic terms of ref. [p!l| . 
That is, we are going to add a^gg of Eq. (|3^) to our Eq. (|25|) and write: 

^W,soft+si ^ ^«,soft ^ -«,si ^ .^1^2, (32) 



where |11| 



aW'^i = -24 In(l-x), 

ag'^^i = -144 ln3(l - x) + ln'^{l - x) (138 - 216 Ijj - 4n/) (33) 
+ ln(l - x) (-132 + 138 - 72 - 4 n/ + 144 C2) • 

Again, one must be very cautious about phenomenological interpretation of this result, 
because of the aforementioned approximate nature of (J^g'^^- 

In Fig.|6| we show the corresponding ii'-factors -K^n°lo^^' -^n°nlIo' ^ proton-proton 
CMS energy of (a) 14 TeV and (b) 2 TeV. By comparing to Fig.^ one can see that the 
formally sub-leading terms clearly dominate the cross section. Figure ^ shows the effect of 
the approximate sub-leading terms on the total cross section. One should bear in mind that 
the agreement of the solid and dotted lines in these two figures is due to the dominance 
of (Tg^g'^^ which has been added both to our ag^g'^°^^ and the approximate ^g^g'^°^^- The 
NNLO corrections shown in Figs. |^ and are quite large and would indicate a very slow 
perturbative convergence of the total cross section. However, it is expected that the use 
of NNLO parton distributions, when they become available, will moderate the size of the 
NNLO terms. 
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j^soft+si^j4 TeV) 




j^soft+sy2 TeV) 



100 120 140 160 180 200 220 240 260 280 300 

My [GeV] 

(a) 




loo 120 140 160 180 200 220 240 260 280 300 

M^^ [GeV] 

(b) 



Figure 6: K factor as defined in Eg. (P0|), including soft and "sub-leading" terms, 
i.e. using Eq. (^) for the dashed and solid curves, and Eq. (^Tj) (with Ug^^ extracted 
from ITiri) for the dotted one. 



a (14 TeV)[pb] 




%0 120 140 160 180 200 220 240 260 280 300 

Mjj [GeV] 




M„ [GeV] 



Figure 7: Cross section a{pp — > H + X) including soft and "sub-leading" terms (c/. 
caption of Fig.|6|). Dash-dotted, dashed, and solid line correspond to LO, NLO, and 
NNLO results obtained from Eq. (p5|). The dotted line is the NNLO result based on 
the approximation given in jl^. (a): = 14 TeV, (b): = 2 TeV. 
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5 Conclusions 



The partonic cross section gg H + X has been calculated to NNLO in the limit of 
soft emission. The resulting hadronic Higgs production rate exhibits well-behaved pertur- 
bative convergence properties. However, the effects of soft emission are expected to be 
significantly smaller than the formally non-leading contributions oc ln*(l — x), for which 
an approximate result has been obtained in ref. |11|. This makes the complete evaluation 
of the NNLO corrections even more imperative. 
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Note added 



The subject of this paper has also been addressed by Catani, de Florian and Grazzini p2[ . 
Their analytical results are in full agreement with ours. Note, however, that their definition 
of the soft limit differs from ours by a global factor of x (cf. Eq.(6) of our paper to Eq.(2) 
of [32|). This means that in |32| some of the sub-leading terms (cf. Eqs.(32),(33)) are 
attributed to the soft contribution, so that the numerical result for our soft limit is not 
directly comparable to their "SV" approximation. Other numerical differences for the final 
results (our agg'''"^^'^^^ , their a^^'-^) are due to the different sets of parton distributions 
and the different treatment of unknown sub- leading terms proportional to ln^(l — x) and 
ln(l — x). 
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